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REPORT 


1. Photoreactor monolith fundamental studies (Y. LUO) 

(a) Monolith reactor models 

Catalytic monoliths have a good prior history in automobile 
exhaust control for emission reduction of carbon monoxide and 
unburned hydrocarbons. Prior reactor modeling of this catalytic 
configuration operating in a thermal (rather than the present 
photochemical ) mode has included either fully developed flow 
patterns or plug flow (flat) velocity patterns. The major exception 
to this claim is the work of Boersma et al (1987) who included the 
more difficult analysis involving the development of the full 
velocity field from the entrance condition of a flat or plug flow 
velocity profile. Long channel monoliths can utilize the fully 
developed flow approximation, since the (nearly) parabolic developed 
flow profile is achieved early in such a channel system. Our 
photocatalytic monolith will use a channel length of six inches 
(commercial auto exhaust configuration) or shorter (to allow for 
more internal illumination) We evaluate this entrance problem 
below. 


(b) Entrance length effects 

The velocity model required depends on the fraction given by the 
ratio of the entrance distance required for parabolic flow 
development divided by the full channel length. A small fraction 
means that almost the entire channel operates under developed flow; 
a ratio closer to unity means the opposite. 

The criteria for estimating the entrance lengths for the velocity 
and concentration fields depend on the Reynolds number (Re) and the 
Schmidt number (Sc). 

For velocity, the entrance length in terms of monolith channel 
diameter, d, is given by 
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Z/d = 0.06 Re 


For concentration fields, the entrance length is given by 
Z/d = 0.06 Re Sc 

For gases, the Schmidt number is unity, thus the entrance lengths 
for development of both concentration and velocity fields is the 
same. 

As the range of Reynolds numbers examined is between 10 and 
150, the velocity entrance length is 0.6 up to 9.0 channel diameters, 
coreesponding to 2% up to 30 % of the channel length. Because the 
illumination field will actually be non-uniform, and most intense 
near the entrance , we expect that a full entrance flow analysis is 
required, since the entrance criteria defined above are from 
literature sources which assume a uniform wall activity, whereas 
the non-uniform photocatalyst wall activity will be strongest in the 
developing flow entrance region. 


(c) Gas flow velocity field 

The velocity flow field is summarized in Table 1. The model 
assumptions include steady state, incompressible flow, Newtonian 
fliuid, axisymmetric velocity profile, and a uniform velocity profile 
upstream of the monolith channel entrance. 

(d) Gas concentration field 

The concentration field equations for the same problem are 
summarized in Table 2. This is a full analysis, including both axial 
and radial dispersion of the contaminant(s) in the gas phase flow. 
The need for axial dispersion, neglected in most previous monolith 
analyses, will be shown in the results section. 

(e) Illumination: uniform vs. non-uniform 

Model calculations were carried out for two circumstances: 
uniform wall intensity and non-uniform intensity, corresponding 
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respectively to uniform and non-uniform wall photocatalyst rate 
constants. 

(f) Results: acetone, 1 -butanol design calculations 

Using the photocataiytic rate constants for acetone and 1 -butanol 
conversion reported earlier in this grant by Peral and Ollis (1992), 
we calculated concentration profiles for acetone and butanol as a 
function of axial position along the monolith (length = 30 channel 
diameters , i e., d = 0.4 cm and L = 30 x 0.4 = 12 cm or about 5 
inches. 

Acetone conversion profiles (Figure 2) indicate a nearly linear 
pattern at all Reynolds numbers (flow rates), consistent with the 

expectation that for a feed concentyration of 712.5 mg/m 3 (SMAC 
level), the reaction kinetics are zero order in acetone (Peral and 
OHis( 1992)). Conversion varied between about 5% per pass at Re = 
150 to about 65% per pass at Re = 10. 

Lower acetone feed concentrations exhibit the curvature expected 
of positive order kinetics (Figure 3; examine curve for feed = 0.1 
SMAC or 71.2 mg/m 3 ) 

Butanol conversion profiles are strongly curved, again as 
expected since the conversion kinetics are nearly first order in 1- 
butanol at the SMAC levell assumed of 121. mg/m3. (Peral and 
Ollis(1 992). 

The same variations of butanol feed concentrations produces the 
same 

C( z )/C 0 curve, reflective of the essentially first order variation of 

rate with concentration at all concnetrations under consideration 
(Figure 4) 

(g) Mass transfer influence 

The computed centerline, mixing cup averaged, and wall 
concentration profile for acetone (Figure 5a) shows very little 
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difference . This result indicates a nearly uniform radial 
concentration profile for acetone, i.e, radial mass transfer rates are 
sufficient and the overalll rate is limited by the photocatalyst wall 
kinetics. 

For the more reactive 1 -butanol, however, the same calculations 
(Fig 5b) indicate that substantial differences of 10-30 % between 
centerline and wall concentrations may exist at any point in the 
channel. This results means the mass transfer resistances are 
appreciable and the presumed requirement for radial mass transfer 
by convective diffusion using a second order radial term is verified 
(Table 2). 

(h) Non-uniform illumination influence 

Calculations assuming an illumination decrease varying inversely 
with distance between the lamp and channel wall are summarized in 
Figures 6a (acetone) and 6b (1 -butanol). A point source lamp is 
positioned at 1, or 10, or 50 channel diameters upstream of the 
monolith (corresponding to 0.4 cm, 4. cm and 20 cm., respectively) 
For acetone, the corresponding conversion per pass for a point 
source lamp changes from 100 % to 35 % to 1 % (Figure 6a). For 
butanol, the corresponding conversions are 100 %, 100 % and 25 % 
(Figure 6b) 

These calculations indicate that for a point source lamp , the 
variation of rate with lamp position is very important. A fuller 
calculation using a more realistic planar light source approximation 
for our lamps is underway and will be reported at a later time in 
this study. 


(i) Conclusions 

The first photocatalytic monolith model has been constructed 
and used to predict concentration profiles for a weakly reactive 
(acetone) and a moderately reactive (butanol) pollutant for various 
conditions of flow (Reynolds number), entrance concentrations, and 
illumination fields. The model appears robust, and we will improve 
it by inclusion of more realistic illumination fields and by 
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confrontation withmeasured illumination and axial concentration 
fields later in this study. 


2. Monolith reactor operation: batch recirculation system (M. Sauer) 


(a) Expermmental set-up 

A recirculating monolith reactor, analogous to that published by 
Suzuki et al of Toyota Laboratories (1991) was constructed of 
transparent plexiglas tubes; a schematic appears in Figure 7. the 
reactor consists of 6 * diameter tubing, two 100 watt near UV 
("black") lights, a variable speed recirculation fan (to allow various 
Reynolds number operations) and a sample injection and 
withdrawal port. 

(b) Results 

(i) Acetone adsorption 

All application of photocatalysis to air purification and 
treatment is likely to involved humidified air. As both acetone and 
water adsorb on the monolith support as well as the photocatalyst 
titanium dioxide phase, and as the support mass is far larger than 
that of the catalyst 10 pm coating, the adsorption isotherms on the 
support need determination. Figure 8 presents the variation of 
acetone mass adsorbed as a function of gas phase concentration of 
acetone at five different injected water contents. (The relative 
humidities were not measurable at this time; an equipment 
modification is underway to measure relative humidity and 
contaminant levels simultaneously) This figure shows the 
importance of relative humidity in influencing the acetone 
partitioning between the gas phase and monolith surface . Our 
reactor models includes this influence explicity. 
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(ii) Acetone kinetics 


Batch experiments in the recirculating apparatus of Figure 8 were 

accomplished at initial concentrations between 100 and 700 mg/m 3 
acetone. Initial rate analysis indicated that the data could be 
fitted to a Langmuir Hinshelwood model, as indicated by the 
linearity of a plot of inverse initial rate vs inverse initial 
concentration (Figure 9) (lowest concentration point ignored, as this 
point was our first run and is not considered particularly reliable) 


(iii) Model development 

The conversion per pass is several percent under appreciable flow 
rates in the recirculation system. Under this circumstance, the 
monolith reactor can be used as a test of kinetic forms over the 
entire course of a batch conversion, since the system as a whole 
behaves as a well stirred reactor. (Later experiments in spring and 
summer of 1993 will explore integral conversion operation at 
slower flow rates) Under this circumstance, we have developeed a 
kinetic model for the reactor system, summarized in Appendix A (pp 
A.1-A.4) The final equation contains only the independent variable , 
time, and the sole dependent variable, C(t). The kinetic model 
includes the following: 

acetone and water adsorption on monolith 

acetone photocatalyst reaction 

system and catalyst volumes 

water photocatalyst inhibition (from Peral et al (1992)) 

constant relative humidity (water of reaction small vs. water for 
humidification) 
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(iv) Comparison: model vs experiment 

The data and model for five initial conditions are summarized in 
Figure lOa-e. The agreement is generally good, and we conclude that 
for acetone conversions, where no appreciable intermediates are 
detectable, the simple model in Appendix A is adequate for reactor 
design purposes. 

Given the adequacy of the model here, it is also useful to consider 
the predicted influence of water vapor variation. Calculations to 
this end are presented in Figure 10. An improved version of this 
model will be developed in spring -sumer of 1993 when we are able 
to determine the water adsorption isotherm experimentally and 
incorporate this more realistic function for the humidification 
influence on both the monolith and the catalyst. 


(v) Dual reactant experiments 

Two runs have been performed to date using dual reactant feed. 
Figure 12a indicates that m-xylene, an aromatic, is converted more 
rapidly than acetonein a two component acetone-xylene feed. This 
result is expected as unsaturates have intermediate to high 
reactivity vs saurated oxygenates. Figure 12b compares the 
acetone vs time concentration profiles for acetone only (triangles) 
and acetone with m-xylene (open squares). This figure indicates 
that a clear competitive inhibition exists when multiple reactants 
are converted simultaneously. This competitive inhibition is found 
to occur also with all other forms of catalysis where a large number 
of reactant molecules compete for a limited number of catalyst 
sites. The next year will include development of two and three 
component contaminant models and confrontation with data such as 
that of Figures 12a and 12b. 
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Table 1 


II. Problem formulation 
Entrance velocity field 


Model assumptions 

(1) . Steady state 

(2) . Incompressible flow 

(3) . Newtonian fluid 

(4) . Axisymmetrical velocity profile 

(5) . Uniform velocity profile before tube entrance 


The dimensionless forms of the governing equations are: 
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Table 1 (cont ? d) 


The only non zero component of the vorticity vector is the 
azimuthal component, which is given by: 


co = 


dV dU 
dZ dr 


The equations we need to solve become: 


r2 az r az & 
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The boundary conditions are: 


Z=0, 0<r<0.5, co=0, \j/=-0.5r 2 
Z>0, r=0, V=co=0, \|/=0 

Z>0, r=0.5, co= 2—, w=-0.125 

3r 2 
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Table 2 


Entrance concentration field 

Surface reaction on the catalyst wall 

The dimensionless governing equation for concentration 
profile in the single channel of monolith reactor is: 

^ + I^ + ^ = Pe^ + PefW 
<3r- 3r dZ? 2 32 2r 3r 


The boundary conditions are: 


3C 
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dr 
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dr 
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El. Reactor analysis for pollutant 
conversions 

• Constant light intensity 

• Point source light intensity 



Fi; 


c/cO 



Computed acetone concentration profiles for Re=10, 50, 100, 150. 

(c0=712.5 mg/m 3 , d=0.4 cm) 


Figure 2a 
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ta — o.icO 
+ O^cO 

* cO 

* 2c0 


Computed profiles for various initial acetone concentrations 
(Re=10, c0=712.5 mg/m 3 , d=0.4 cm) 


Figure 2b 
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c/cO 
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Computed profiles for various initial 1-butanol concentrations 
(Re=10, c0=121.0 mg/m 3 , d=0.4 cm) 


Figure 4 
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z(l/d) 

Computed acetone wall, centerline, and average concentrations 
(Re=10, c0=712.5 mg/m 3 , d=0.4 cm) 

Figure 5a 
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Computed 1-butanol wall, centerline, and average concentrations 
(Re=10, c0=121.0 mg/m 3 , d=0.4 cm) 

Figure 5b 
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Acetone concentration profiles 

(c0=712.5 mg/m 3 , d=0.4cm, Osram HBO 500W point source lamp) 


Figure 6a 
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1-Butanol concentration profiles 

(c0=121.0 mg/m 3 , d=0.4cm, Osram HBO 500W point source lamp) 


Figure 6b 



Experimental monolith reactor system setup 



Figure 7 
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M_a(mg) 


Acetone adsorption 

• Mass of acetone adsorbed on monolith (M A ) verses 
steady state acetone gas phase concentration (C A ) 


M_a vs Ca 



Figure 8 





Model/Experiment comparison 


Model prediction for Co — lOO^f experimental data 


Cg vs Time 

Cp = 100 mg/m * 
212 mg water 



. Model 


fs *1 os Os 0»1 *0 *Oi S 3 S 3 ts ts 


time( min) 


Figure 10a 

Model prediction for Cq = 250^? experimental data 


Cg vs Time 
Co™ 250 mg/ nr* 
212 mg water 



. Model 
Exp . Data 
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Figure 10b 


Cg(mglm?) ® Cg(mg!n?) 


Model prediction for Cq = 400^f experimental data 



time( min ) 

Figure 10c 

Model prediction for Cq = 700^ experimental data 

7T2> 



time(min) 


Figure lOd 
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• Model prediction for Cq — lOOO^f experimental data 

* / L 



time( min) 

Figure lOe 
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• Short time model prediction of water inhibition effect 



T 


t24 mg water 
MS mg water 
M2 mg water 
\Q6 mg water 
i.2 mg water 


Figure 11a 


Long time model prediction of water inhibition effect 


Cg vs Time 
Co ss 400 mg/m ^ 



424 mg water 
31S mg water 
212 mg water 
106 mg water 
4J2 mg water 


Figure lib 
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Two reactant experiment 

• Acetone/m-Xylene mixed feed experimental result 


Concentration vs Time 
Acetone/n 




Im-Xylene 

= 400 mg/m 3 
mg water 



d Acetone 
o— m-Xylene 


Figure 12a 
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Acetone Concentration (mg/m 


• Comparison of Acetone Single Feed and Mixed Feed 
Runs 


Concentration vs Time 
Acetone 

Co = 400 mg/m * 

212 mg water 



Time(min) 


□ Mixed Feed 

a Single Reactant 


Figure 12b 


28 



Appendix A 


Model development 


ft A reactant species mass balance assuming the recycle loop is plug flow (PF), 


Assumptions (PF) : 

• Plug flow 

• No reaction occuring 

• Isothermal 


If we assume a well mixed recycle loop, (i. e. C g ^ /(z)), with A a L 

= 9(02 _ Cl) 

Equation 2 is the equation for a CST, with C $ = C\ 

■Jb A reactant species mass balance on the CSTR gives, 

fr r , kKC R V^ _ dC R dM R 

9(Cl ~ C2) " TTxcT * dt + dt 

Assumptions (CSTR): 


= V„ 


( 2 ) 


( 3 ) 


• Neglect mass transfer effects 

• Reaction follows a simple LH rate form 

• Well mixed system 

• C R « = C 

• Isothermal 


A. 1 
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Combining Equations 2 and 3, 

kKCuV?* 

1 + KC r 




dC a \ cIMr 


+ 


dt 


(4) 


V. + V R = V s (5) 

The conversion per pass for Cq = lOO^f is ~ 0.45%. For a low per pass conversion 


assume (Cr ~ C*) = C g . This assumption also gives, 


3 dt 

Combining Equations 4 and 6, 




dt 


dt 


( 6 ) 


kKC g V c aci 
1 + KC 9 


g v c = y dC 9 ■ dM R 


'9 


dt 


+ 


dt 


(7) 


Equation 7 has no fit parameters. Substituting the Langmuir isotherm for ace- 
tone, 

« W >» t I / IX # * \ 

( 8 ) 


kKC g V c act _ y dCg L d_ ( i± a K a C, 


1 + 3 dt + dt \1 + + K w C t 

The differentiation £ ( j+rf^t c z ) gives ( assumin S C ™ ~ constant ), 

HaK A^g 


dt \1 + K^Cg + K W C W 

I I^aKa 

\1 + K A C g + K W C W 

\l A K\Cg \ dCg 

(1 +K A C g + K w C w y) dt 

Substituting Equation 9 into Equation 8, 


(9) 


kKCgV^ 


1 + KC q 


— I y a + 


HaKa 


t^A-^A.^9 


dC n 


1 +K A C g + K w C w (l + K A C g + K w C w ) 2 J dt 


( 10 ) 


A .2 
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Rearranging Equation 10 and integrating, 
f kKV^dt = 

- C ( i + *)( v ' + 1 

PaK\C b 


PaK a 


+ K A C g + K W C 

da 


W 


(1 + K A C g + K W C W )*J 3 

Evaluation of tlie integrals (assuming C w = constant ), 
kKV^t = 


(ii) 


( 

-(= 


PaKa 
1 + K W C V 


In 


' CM 

c g (t ) 


+ KV„(C g ( 0) - c,(0) 


PaKa \ ]n (l + K A C g (0) + K„C. 


w 


,1 + K W C 

l^A 


1 + K A C g (t) + K W C 
1 


w 


,1 + K A C,(0) + K W C„ 1 + K A C g (t) + K m C w 


(12) 


Separating the water inhibition effect from the reaction rate constant ( k ) leads 
to a reaction rate constant (fc 0 ) that is independent of water effects. 


k = 


1 + k wa ci 


(13) 


The water inhibition expression for acetone and the values for K Wa and a were 
obtained from Peral and Ollis, and were assumed to be applicable in the current 
study. The water inhibition effect will be examined with the monolith system in 
future experiments. Substituting values for k,K WA and a into Equation 13 gives, 


h = 0.77 


mg 


cm 3 cat — min 


(14) 


Where, 


K Wa = 9.6 x 10~ 7 


m 


mg 


a = 1.674 


C m = 11509 


mg 


m* 


A. 3 
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The resulting model equation is, 


fco KV c act 


(' 
-G 


HaK A 


1 + K W C W 
VaK a 


'CM' 


d+K w C, 
—HaKK w C, 


M# +jrW) " CjW) 

\ Iti ( l±KACMj L K^ 

’ J \ 1 + K A C g {t) + K W C W 


K 1 + KaCM + K w c w 1 + K a CM 


Equation 15 utilizes no fit parameters: 

• Acetone adsorption Langmuir fit ==$■ (i A ,K A ,K w 

• Acetone kinetic parameters ==> k,K 

• System physical constants => V act V a 

• Dr. Peral’s water inhibition constants => K Wjl ,a 

• Water independent reaction rate constant =~> jk 0 

• Assume C w = constant, C ff (0) = constant (initial gas phase acetone 
tration) 


A. 4 


(15) 


concen- 
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